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Localized astronomical sources like a double stellar system, rotating neutron star, or a massive 
black hole at the center of the Milky Way emit periodic gravitational waves. For a long time 
only a far-zone contribution of gravitational fields of the localized sources (plane-wave-front 
approximation) were a matter of theoretical analysis. We demonstrate how this analysis can 
be extended to take into account near-zone and intermediate-zone contributions as well. The 
formalism is used to calculate gravitational-wave corrections to the Shapiro time delay in 
binary pulsars and low-frequency (LF) pulsar timing noise produced by an ensemble of double 
stars in our galaxy. 



1 Introduction 

The extremely high precision of current pulsar timing observations demands a better theoretical 
treatment of secondary effects in the propagation of electromagnetic signals in variable gravita- 
tional fields of oscillating and precessing stars, or binary systems. Especially important is the 
problem of propagation of light rays in the field of gravitational waves emitted by a localized 
source of gravitational radiation. A consistent approach for a complete and exhaustive solution 
of this problem has been developed in papers ^ra in the first post-Minkowskian approximation 
of General Relativity. We have demonstrated for the first time that the equations of light prop- 
agation in the retarded gravitational field of an arbitrary localized source emitting gravitational 
waves can be integrated exactly in closed form. It allows to examine the influence of the time- 
dependent gravitational field on the light propagation not only in the wave zone but also in 
cases when light passes through the intermediate and near zones of the source. We have ob- 
tained explicit analytic expressions for light deflection and integrated time delay (Shapiro effect) 
accounting for all possible retardation effects and arbitrary relative locations of the source of 
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gravitational waves, the source of light rays , and the observer. Coordinate dependent terms in 
the expressions for observable quantities were singled out and used for the unique interpretation 
of observable effects. Our exploration essentially extends previous results regarding propagation 
of light rays in the field of a plane monochromatic gravitational wave (see, e.g., papers 
and significantly surpasses theoretical approaches of other authors B~t3. In the present paper 
we briefly discuss the developed formalism and apply it to describe the Shapiro time delay in 
binary pulsars as well as to estimate the intensity of LF timing noise produced by an ensemble 
of double stars of our galaxy. 

2 Propagation of light in time-dependent gravitational fields 

We solve the Einstein equations in the first post-Minkowskian approximation where the metric 
tensor g a p is decomposed linearly into the Minkowski metric r) a /3 and a small perturbation /i Q/ g, 

9a(3 = Va/3 + h aj3 . (1) 

The Einstein equations for h a p in harmonic gauge read|| 

□ Kpit, x) = -16nS a p(t, x) , where S a p(t, x) = T a/3 (t, x) - ^q af3 T x x (t, x) . (2) 
Herein T a/3 is the tensor of energy- momentum of a system of massive particles 111 

N 

T^(i,x) = £ ff(t)8 (x - x«(t)) , ff(t) = m al ~\t)u a a it)u^t) , (3) 

0=1 

where t is coordinate time, x = x l = (x 1 ,^ 2 ,^ 3 ) denotes spatial coordinates of a current point 
in space, m a and x a (i) are the rest mass and spatial coordinates of the a-th particle, v a (i) = 
dx a (t)/dt, J a (t) = [1 — ^a(^)] -1 ^ 2 1S the Lorentz factor, u%(t) = {j a (t), 7a(*)v a (i)} is the four- 
velocity of the a-th particle, 5(x) is the 3-dimensional Dirac delta-function. The solution of @ 
is the retarded Lienard- Wiechert tensor potential 
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^(t,x) = 4 V v; ; 2 ' aX ;\ J , (4) 



a=l 



l-a(s) 



where the retarded time s = s(i,x) for the a-th body is a solution of the light-cone equation 
s + |x — x a (s)| = t. Here it is assumed that the field is measured at time t and at the point x. 

We consider the motion of a light ray (photon) in the background gravitational field described 
by the metric (|j) . The motion of the photon is defined by solving equations of a light eeodesic 
(see Fig |l|). The original equations of propagation of light rays are rather complicated □. They 
can be simplified and reduced to a form which allows to resort to a special approximation 
method for their integration. The reduced equations and details of the integration procedure 
can be found in our paper H. In what follows we discuss some of physical applications of the 
mathematical technique. 



3 Shapiro time delay in binary pulsars 

Damour & Esposito-Farese El argued about the contribution of gravitational waves emitted by 
a binary pulsar to the observed Shapiro time delay. As a consequence of the main result of their 
work they claimed that the contribution probably should be small. However, they were not able 
to present conclusive mathematical arguments which are given in this section. 

^In what follows we use geometrical units in which G = c = 1. 
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Figure 1: Illustration of the light-ray's propagation history. The light ray is emitted at the instant of time to at 
the point xo and arrives at the point of observation x at the instant of time t. Light-deflecting bodies moves along 
accelerated world lines during the time of propagation of the light ray; their velocities at some intermediate instant 
of time are shown by black arrows. In the absence of the light-ray-deflecting bodies the light ray would propagate 
along an unperturbed path (dashed line) which is a straight line passing through the points of emission, x U; and 
observation, x. Direction of the unperturbed path is determined by the unit vector K = — (x — xo)/|x — xo|. 
In the presence of the light-ray-deflecting bodies the light ray propagates along the perturbed path (solid line). 
The perturbed trajectory of the light ray is bent and twisted due to the gravito-electric (mass-induced) and 
gravito-magnetic (velocity-induced) fields of the bodies. The initial boundary condition for the equation of light 
propagation is determined by the unit vector k defined at past null infinity by means of a dynamical backward- 
in-time prolongation (dotted line) of the perturbed trajectory of light from the point of emission xo in such a way 
that the tangent vector of the prolongated trajectory coincides with that of the perturbed light-ray's trajectory 
at the point of emission. The relationship between unit vectors k and K includes relativistic bending of light. 



The general formalism of our paper a yields for the time of propagation of the light-ray 

t-t = |x-x | +A(t,t ) , (5) 

where |x — xo| is the Euclidean distance between the points of emission, xo, and observation, 
x, of the photon, and A(t, to) is the Shapiro time delay due to the gravitational field of moving 
bodies 
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A(Mo) = -k a kpB a P(T)--k a kpB a P(T Q ) = 2j2m a B a (s,s ) , (6) 



r [l-k-v a (C)] 2 d( 

where the retarded time s is obtained by solving the equation of light cone for the time of 
observation, the time sq is found by solving the same equation written down for the time of 
emission of photon so + |xo — x a (so)| = to, and t* is the (constant) time of the closest approach 
of the light ray to the origin of the coordinate system used for calculations. Choosing the origin 
of the coordinate system to be the barycenter of the binary pulsar, expanding (0) in powers of 
v a /c, and performing integration result ini 

A(Mo) = -2m c (l-k-v c )ln[(r + k-r)-v c T+(k-r)(k-v c )] , (8) 

where m c is the mass of the pulsar's companion, r is the radial distance between the pulsar 
and companion, v c is the companion's barycentric velocity, and all quantities are taken at the 
barycentric time of emission of the pulsar's radio pulse. Formula (@) apparently demonstrates 
that the velocity-dependent terms present in the Shapiro effect are small and only marginally 
detectable in timing of binary pulsars. This proves the intuitive guess of the paperEJ. Formula S) 



can be implemented in the parameterized post-Keplerian formalism0lll. Its explicit formulation 
is given ini. In case of circular orbits it reads 

2Gm, 
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(9) 



where m p is the pulsar's mass, ip is the orbital phase, ipo is constant, x is the projected semi- 
major axis of the orbit, Pb is the orbital period, and we have restored G and the speed of light 
c for convenience. The correction to the Shapiro delay given in @ may be observed in binary 
pulsars with orbits visible nearly edgewise. 



4 Low-frequency pulsar timing noise 

Each double star in our galaxy emits rather weak LF gravitational waves. However, an ensemble 
of double stars may produce a noticeable amount of the gravitational radiation which could be 
detected through the analysis of the long-term behavior of the timing residuals (see Fig. 2). 
Time variations caused by quadrupolar gravitational waves from a double star read asi 
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(10) 

(11) 



where k l is the unit vector from the pulsar towards the observer, r is the distance 'observer- 
double star'; ro is the distance 'pulsar-double star' (ro = (r 2 + R 2 — 2ri?cos 9) 1 / 2 ); 8 is the angle 
between the directions 'observer-double star' and 'observer-pulsar'; 9q is the angle between the 
directions 'pulsar-double star' and 'observer-pulsar'; Xf^ is the transverse-traceless quadrupole 
moment of the double star; A(9,t — r) is a regular function. If one takes a particular event of a 
circular orbit the formula (|ll]) is brought to the form 
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where cp = tot + a is the orbital phase of the double star, a is constant, to = 2tiv is the orbital 
frequency, O is the longitude of the ascending node, and M. = (mim2) 3 / 5 (mi + ma)" 1 ' 5 is the 
"chirp" mass of the double star. 

An ensemble of double stars generates the gravitational wave noise with the autocovariance 



function H(t) defined as the ensemble average of products of the time delay (10) 

K(t) = < A(t,t-R) A(t + r,t- R + t) > , 



(13) 



where we have used an approximation to = t — R. The ensemble averaging is defined with respect 
to the set of random variables 0, i, ipo, v, 9, and r j]. We assume that the variables O, i, ipo are 
distributed uniformly. The probability distribution function of other variables of the ensemble 
is determined by main-sequence galactic double stars which occupy the frequency bandll3 from 
10~ 13 to 10 -5 Hz with the frequency and spatial distributions F(y) and n(r) respectively f| 
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c We note that variables ro and do are not independent and can be expressed through r and 8 
d Function F{u) is nomalized to unity, and n(r) = n{6, r) is normalized to the total number of double stars in 
the ensemble. 
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Figure 2: A model of an ensemble of double systems producing LF gravitational noise. The minimal distance 
observer-double star is denoted as r m in and the minimal angle between directions to the pulsar and a double 
star is denoted as 8 m i„. The distance from the observer to the pulsar is R. Concentration, n, of double stars is 
assumed to be constant. Then, r mi „ ~ n~ 1//3 , 9 min ~ r min /R, and n(r) = 3n/(4-7ri? 3 ). 

where it has been assumed that F(y) and n(r) are statistically independent. Function k(i>) is 
almost constant over a wide range of frequencies. Hence, we take F(v) ~ 0.054/za One can 
prove that the double-star-gravitational-wave noise is stationary with spectrum S(u) defined by 
the Fourier transform of the autocovariance function. The spectrum has the following specific 
properties 

• the noise consists of seven components having spectra S(v) ~ i/~ 5 / 3 ~ a (a = 4- 6); 

• the amplitude of any component of the noise crucially depends on the spatial distribution 
of galactic double stars between the observer and the pulsar as well as in the neighborhood 
of each (e.g. a pulsar in a globular cluster); 

• the component of S(v) with a = is produced by the far zone gravitational field (~ r _1 ) 
of double stars (pure gravitational waves) and is regular everywhere including the line of 
sight of observer towards the pulsar, when 9 = 0; 

• the component of S(u) with a = 2 is due to the semi-far zone gravitational field (~ r~ 2 ) 
of double stars and diverges along the line of sight of observer towards the pulsar fl; 

• the components of S(u) with a = 4, 6 are due to the intermediate zone (~ r~ 3 ) and near 
zone (~ r -4 ) gravitational fields of double stars and are regular everywhere. 

In order to give an estimate of the expected amount of the gravitational wave noise in pulsar 
timing we use a simple model of uniform distribution of double stars in the spherical domain 
between the pulsar and observer at the Earth (see Fig. 2 for more details). The strength of the 



e This divergency of the spectrum is cJpsely related to the effect of augmentation of the time delay by gravita- 
tional lensing by the intervening matter □ 



timing noise is most conveniently characterized by the, so-called, a z statistic 113 defined as 



'.(t) = < cl(r) >V2 , (15) 

where 03(7") is a random variable being proportional to the statistical fluctuations of the third 
time derivative of the pulsar's timing residuals, r(t), and r is the total span of observational 
time. In the case of a single double star, r(t) = A(t,t — R). In what follows we simplify the 
problem by accounting for only the line-of-sight divergent component of A(i, t — R) given in the 
second line of (|l2|). Performing integration over the ensemble's variables and over the sperical 
volume {] shown in Fig. 2 gives 

rl/r 

< c|(r) >~ (2tt) 6 / u s S(v)dv , with S(u) ~ z^ 11/3 , so that (16) 



(0 /3 fe)(r^)fe) 1/3 - <") 

This limit is lower than that presently accessible in the timing of millisecond pulsars. However, 
we emphasize that the estimate (|i~7|): 1) is model dependent, and 2) contribution of several other 
components of A(t, t — R) were neglected. It is highly desirable to repeat our calculations with 
more realistic ensemble of double star distributions and the complete account for all constituents 
of A(t,t — R). We argue that proceedingin this way the origin of the incomprehensible 'red' 
timing noise discovered in PSR B1937+21E3 may be explained as a result of the LF gravitational 
wave noise from double stars in our galaxy - a challenge both for theorists and observers. 
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